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Abstract 

We come back to a non linear integral equation satisfied by the function H, which is 

■ distinct from the classical ii-equation. Established for the first time by Busbridge 
(1955), it appeared occasionally in the literature since then. First of all, this equation 
is generalized over the whole complex plane using the method of residues. Then its 
counterpart in a finite slab is derived; it consists in two series of integral equations 
for the X- and y-functions. These integral equations are finally applied to the 
solution of the albedo problem in a slab. 
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C3 ' 1 Introduction 



This paper originates in a work by Busbridge [1] on the calculation of the 
classical X- and y-functions of transfer theory. A non linear integral equation 
for the if -function is proved in this article, that is not the classical if -equation: 
see Eqs. (5.6)-(5.7) in [1], hereafter Eqs. (27) and (29). 

In the second of a series of papers submitted for publication [2,3,4] , we deduced 
Busbridge's equation from the properties of the iV-function, which can be seen 
as the "heart" of the ii-function. The simplifying role of Busbridge's equation 
for solving the singular if-equation was stressed in Sec. Ill- A of [3]. It also 
clarifies the solution of the albedo problem, yielding the surface field from the 
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internal field of a half-space illuminated from outside. It thus provides the 
link between the internal and the surface solution to the albedo problem in a 
half-space. 

Considering the importance of Busbridge's equation in semi-infinite geometry, 
it is tempting to generalize it in the finite case and examine the consequence 
it has on the albedo problem in a finite slab: this is the object of the present 
article. After a reminder on the classical functions H, X and Y (Sec. 2), we 
come back in Sec. 3 to Busbridge's i?-equation. This equation is generalized 
to the complex plane and proved by the method of residues. This approach is 
then continued in a finite slab, which leads to two series of integral equations 
for the X- and Y- functions (Sec. 4). An application to the slab albedo problem 
is proposed in Sec. 5. 



2 Classical integral equations for the H-, X- and F-functions 

We consider a slab of constant volumic albedo a e]0, 1[ and optical thickness 
b (0 < b < +oo), where light scattering is isotropic and monochromatic. The 
auxiliary functions describing the radiation field depend on the a-parameter 
in a half-space (b = +oo), and on the a- and 6-parameters in a finite slab. 
This dependence is always explicitly mentioned. 

2.1 Integral equations for the H -function in a half-space 

We introduce the dispersion function for the adopted scattering law 
~w a, f +1 du 

T(a,z) = l + -z (z?±l), (1) 

Z J-i u — z 

which describes the multiple scattering of photons in an unbounded medium. 
This function is defined in the whole complex plane C, except at ±1, provided 
that the integral is calculated in the sense of the Cauchy principal value when 
z e] — 1,+1[. When < a < 1 the characteristic equation T(a,z) = has 
four non zero roots in C, namely two pairs of opposite real numbers since 
the T-function is even [5,6]. They are denoted by ±l/k (0 < k < 1) and 
±1/K (0 < K < 1). The important root k in ]0, 1[ is calculated by solving the 
transcendental equation 

T(a,l/*) = l-~ln[±! = (0<*<1), (2) 
which has an exact solution k — k(a) [7]. 

The if-function may be defined as the unique solution, analytic in the right 
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complex half-plane, of any of the following integral equations 

H(a,z)[l-Ua,z)] = l (z^-1), (3) 

T(a, z)H(a, z) = 1 - £„(a, -z) (z ^ ±1), (4) 

[l-£ H (a,z)][l-Ua,-z)]=T(a,z) (z^±l), (5) 



where 



n c\ii 
Ua,z) = -z H(a,u)—— (z ? -1) (6) 
2 Jo u + z 

(the integral is a Cauchy principal value for z e] — 1, 0[). 

Putting 2 = +1/A; and z = +1/K into Eq. (4) and observing that the H- 

function does not diverge at these points, one obtains 1 — £ H (a, —1/k) = 

1 — £ w (a, —1/K) = 0. Equation (3) then shows that H diverges at z — — 1/k 

and z — —1/K. This equation also defines the extension of the if -function in 

C \ {-1/Jfe, -1, -1/K}, and we have for any 2 iiiC \ {-1/k, -1, -1/-K"} 

— , = l-^ H {a,z) = \--z H(a,u)——. 7 

H(a,z) 2 Jo u + z 

Equation (4) reads also, for any z G C \ {— ±1, —1/K} 

T{a,z)H{a,z) = l + -z H(a,u)——, (8) 

2 JO M — z 

which allows to extend H by continuity at z — —1: H(a, —1) = 0. Finally, it 
follows from Eqs. (3) and (5) that for any zinC \ {±l/k, ±1, ±1/K} 

T(a, z)H(a, z)H(a, -z) = 1. (9) 

Equations (7)- (9) are the regular if -equation [8], the singular if-equation [9], 
and the factorization relation [10] respectively. We refer here to those papers 
where, to our knowledge, each equation can be found for the first time. 
The behavior of the if -function at infinity is required in Sec. 3. We recall that 

if(a,oo) = — (10) 
V 1 — a 

as can be seen by putting z — > 00 in Eqs. (4) and (5), observing that T(a, 00) = 
1 - a. 



2.2 Integral equations for the functions X and Y in a finite slab 

The X- and F-functions are defined in C* = C \ {0} and related by 
X(a, b, z) = Y(a, b, — z) exp(— b/z), Y(a,b, z) — X(a,b, —z) exp(— 6/ z). {11) 
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They are the unique solution, analytic in the right complex half-plane, of 
coupled integral equations involving the auxiliary functions 



£ x (a,b,z) = -z X(a,b,u)——, (12) 
2 Jo u+ z 

£ Y (a,b,z) = °-z ( Y(a,b,u) — ^— , (13) 

A Jo U + Z 

defined over C \ { — 1}. These equations are either the regular X- and In- 
equations [11] 

X(a, 6, z)[l - £ x (a, b, z)\ + Y(a, b, z)£ Y (a, b,z) = l (z ^ 0, -1), (14) 



X(a, b, z)£ Y (a, b, -z) + Y(a, b, z)[l - £ x (a, b, —z)] = exp(-b/z) (z ^ 0, +1), 

(15) 

or the singular X- and F-equations [9] 
T(a, z)X(a, b, z) = 1 — £ x (a, b, — z) — £ y (a, b, z) exp(— b/z) (z^0,±l), (16) 

T(a, z)Y(a, b, z) = [l-^(a, b, z)\ exp(-&/z)-£ y (a, b, -z) (z ^ 0, ±1). (17) 
We have furthermore the relation [12] 

[1 - £ x (a,b,z)][l - €x(a,b, -z)\ - £ Y (a, b, z)£ Y (a, b, -z) = T(a,z) (z ^ ±1), 

(18) 

which generalizes the factorization relation (9) for semi-infinite media. 
Let b go to infinity and suppose that Re(z) > 0, where Ke(z) is the real part 
of the complex number z. Then X(a,b,z) — > H(a,z) and Y(a, b, z) — > 0, so 
that Eqs. (14), (16) and (18) yield Eqs. (3)-(5) again respectively. Equations 
(15) and (17) reduce to = 0. 

The behavior of the functions X and Y in a neighborhood of and at in- 
finity follows immediately from Eqs. (16) and (17). Let z tends to in these 
equations. Then T(a,z) and 1 — £ x (a,b,z) tend to 1, and £ Y (a, b, z) tends to 
(a/2)(3_ 1 (a,b)z, where /3-i(a, b) is the moment of order -1 of the F-function, 
viz. 

/3_i(a,6)= / Y(a,b,u)— . (19) 
Jo u 

Consequently 

X(a,b,z) ~ l-^f3_ 1 ( a ,b)zexp(-b/z) (z -> 0), (20) 
Y(a,b,z) ~ exp(-b/z) + ^(3_ 1 (a,b)z (*->()). (21) 
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These relations show that X(a,b,z) — > oo like — (a/2)/3_i(a, b)zexp(— b/z) 
and y(a, 6, z) — > oo like exp(— ft/z) when z — ► in the left complex half-plane 
[Re(z) < 0], and X(a,b,z) — > 1 and F(a, b, z) — > when z — >■ in the right 
complex half-plane [Re(z) > 0]. 
Now letting z — > oo into Eqs. (12)-(13) yields 

£ x (a, 6, oo) = ^a {a, b), £ Y (a, b, oo) = |/3o(a, b), (22) 

where a (a,b) and (3 (a,b) are the moments of order of the X- and F- 
functions 

a (a,b) = / X(a, 6, w) du, (3o(a, b) = / Y(a,b,u)du. (23) 
Jo Jo 



3 Busbridge's integral equation for the //-function 



This integral equation was established in 1955 by Busbridge [1] as an auxiliary 
result for the calculation of the X- and F-equations. This first appearance 
was so discreet that Busbridge didn't include her lemma in [5]. Busbridge's 
equation was "rediscovered" by Andreasian and Danielian [13], and Das [14] 
wrote a version which needs simplification. It appears also in the Appendix D 
of [15]. The domain of this equation has been extended to the real axis in [16] 
and to the complex plane in [17] and [3]. In the two references mentioned at 
last, it is noted that Busbridge's equation clarifies the albedo problem, since 
it allows deriving the surface source function from the internal one, a result 
we intend to generalize in a finite slab (Sec. 5). 

Here we summarize the method used in [16] and [17], which is based on the 
residue theorem. It is well-suited to the main aim of the present paper - the 
generalization of Busbridge's equation in a finite slab - since it still holds in 
the finite case. It consists in integrating the function 

f:z'eC^H{aMz')—^— i 

1 — zz' 

along the contour of Fig. Al in [16], replacing H(a, +1/V) in the left complex 
half-plane by l/[T(a, +l/z')H(a, —l/z')\. Results involve the function 

9(a, v) = — — ^- — (0 < v < 1), (24) 

T z (a, v) + yKav/2y 

which arises when skirting the ] — oo, — 1[ cut, and the coefficients 
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S(a,k) = 1 



ak 2 



(26) 



(k 2 + a - 1)(1 - k 2 y 
related to the pole z = —k within the contour. We recall that k = k(a) is the 
unique solution in ]0, 1[ of the transcendental Eq. (2). 
It follows from the residue theorem that, for z G C \ {[—1, 0[U{ — l/k}}, 



H(a,z) = 1 + 



R(a, k) 



kz 



H(a,+l/k)l + kz 



+ 



°-z [\g/H)(a,v)-^- (27) 
2 Jo v + z 



where (g/H)(a,v) = g(a,v)/H(a,v) is a notation we adopt hereafter. 
This equation is still valid at z — — 1, provided that H(a, z) is replaced by its 
extension by continuity at z = —1, which is 0. Hence the exact value of the 
integral 

o Z* 1 / / TT w , dv „ R(a,k) k , . 

3 / (,///)(»,„)_ = j.-^-L^— (28) 

For z = u e] — 1,0[, the left-hand side of Eq. (27) must be replaced by 
(gT/H)(a, —u) because of the presence of the pole 1/m on the cut ] — oo, — 1[. 
Moreover, the integral on the right-hand side of Eq. (27) has to be calculated 
in the sense of the Cauchy principal value (symbol /). Hence 



(gT/H)(a,-u) = l + 



R(a, k) 



ku 



a 

j u 

H(a,+l/k)l + ku 2 



-f\g/H)(a,v)-^-, (29) 
Jo V + u 



which yields Eq. (28) again as u — > — 1 on the right. 

Finally, Eq. (27) can be extended by continuity at z = — l/k by 



\g/H)(a,v)- d ' 
o 1 



kv 



-1- 



R(a, k) 



H(a,+l/k) 



S(a, k) + 



1 H'(a,+l/k) 
kH(a,+l/k) 



, (30) 



where H'(a,+l/k) is the derivative of the if -function at +l/k. It can be 
deduced from the derivative of the if-equation (7) with respect to z. 
Letting z — > oo into Eq. (27) and using Eq. (10), we obtain 



y\g/H)(a,v)dv 



1 



R(a, k) 
H(a,+l/kY 



(31) 



A new series of integral equations can be deduced from the previous one by 
substracting the above equation to Eqs. (27)-(30): see Sec. III-B of [3]. 
The equations demonstrated by Busbridge [1] are (27) and (29). Equation (27) 
is also in Refs. [13,14,15,16,17]. Equation (29) was retrieved in [16] and [17], 
which contain also Eqs. (28), (30) and (31). The whole set of equations and 
their immediate consequences is derived in [3] from a preliminary study of the 
iV-function [2]. 

Our focus now is to write the counterparts of Eqs. (27)-(31) in a finite slab 
(Sec. 4), then to investigate some application for the solution of the albedo 
problem (Sec. 5). 
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4 The finite case: integral equations for the X- and F-functions 



These equations can be proved by the method of residues, using the same con- 
tour as in the semi-infinite case (Fig. Al in [16]). Just integrate the functions 

z' e C^X(a,b,l/z') ' 



1-zz' 



z' e C -> Y(a, b, l/z') , 1 , 

1 — zz' 

along this contour, writing in the left complex half-plane 



-£„(„, 6, (32) 
1 + zz J 



/N1 exp(6z / ) 

+ 22/ 



^(a,6,-l//)-^| dz', (33) 
1 — 22 J 



due to Eqs. (16)-(17). This transformation overcomes the difficulty arising 
from the divergence of the functions z 1 — > X(a, b, l/z') and z' — > F(a, 6, 1/V) 
when z' — > oo in the left complex half-plane [see Eqs. (20)-(21)]. One obtains, 
forzeC\{[-l,+l]U{±l/ib}} 

X(a, 6, z) + — -f y (a, 6, z) exp(-b z) = — 

1 {a, z) 1 (a, z) 

I kz kz ii 

= 1 + R(a, k)\[l- £ x (a, b, +l/k)\- — - £ Y (a, b, +l/k) exp(-kb)- 



1 • /.•; " ' 1 — kz 



+ °-z f 1 g(a, v)\[l- Ua, b, v)]-±- - £ y (a, 6, v f M bM ) dv, (34) 
2 Jo [ v + z v — z J 



Y(a, b, z) - -^—[1 - ^(a, b, z)} exp(-6/z) = - U °> b > " ] 
1 (a, z) 1 (a, z) 

= R(a, k) |[1 - £ x (a, 6, +l/fc)] exp(-A;6) T ^- - £ y (a, 6, +1/*;)^-^- j 

+ ^ /^(a,*) {[1 - Ua,b,v f M ~ b/v) - Z Y ( a ,b,v)^—\ dv. (35) 
2 Jo [ v — z v + z J 

At 2 = ±1, Eqs. (34)-(35) still hold, provided that 1/T(a,z) is replaced by 
in both members of the first equality. This simplification has to be done in 
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the left-hand side for z — +1, and in the right-hand side for z = — 1. One has, 
at z — +1 



X(a,b,+1) = 1 

+ R(a, k) |[1 - k(a, 6, + V*)] j"^ ~ 6 ' ex P(-^)^} 

+ £ C 9(a, v)\[l- U«, b, v)]^— - &a, b, v fM^M \ dw> (36 ) 
2 jo [ v + 1 f — 1 J 

y(a, 6, +1) = i?(a, k) |[1 - £ x (a, 6, +1/A;)] exp(-fc&)^-£ 

-^(a,6,+l/A;): 



1 + fc_ 

+■ ? f 9(a, v)\[l- Ua, b, v) f M ~ b l v) - £ y (a, b,v)^—\ dv, (37) 
2 Jo [ v — 1 f+lj 



and at z = — 1 



= 1- i?(a, k) ( [1 - £ x (a, 6, +1/A;)]— ^— - £ y (a, 6, +l/k) exp(-kb) 



l-k " v ' ' ' ' cy '1 + JfcJ 

(38) 



a - f g(a, v)\[l- Ua, b, v)}^— - £ y (a, b, v fM b/v) ) 
2 Jo [ v — 1 f + lj 



= R(a, k) ( [1 - £ x (a, 6, +!/£;)] exp(-A;&) — ^- - £ y (a, 6, +1/A;)- 



1 + ife ^ v ' ' ' '1-ife 



+ ? T ff(a, v)\[l- Ua, b, v) f M ^ v) - £ y (a, b,v)^-\ dv. 
2 Jo [ v + 1 f — 1 J 

If z = u e] - 1, 0[U]0, +1[, Eqs. (34)-(35) must be replaced by 



(39) 



Y(w)X(a,6,w) = 1 

+ R(a, k) |[1 - ^x(a, 6, +1/*)]^-^- - £ y (a, 6, exp^^^- j 

- (^T)(a, «) {Y(-«) [1 - ^(a, 6, -«)] + Y(«)^(a, 6, «) exp(-6/ M )} 

+ « A(a,,) f[l -M^M)^ 1 ^! dt;, (40) 

2 Jo I v + u v — u J 
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Y(u)Y(a,b,u) = 
R(a, k)\[l- £ x (a, b, +l/k)\ exp(-kb)- £ y (a, 6, +l/k) 



1 — ku 1 + ku J 

+ (gT)(a,u) {Y(u)[l - Z x (a,b,u)]exp(-b/u) + Y(- U )£ Y (a,b, -«)} 

+ v)\[l- U", b, v) f M ~ bM - £ y (a, M)~M dt;, (41) 

2 Jo [ f — u v + u J 

where Y denotes the Heaviside function: Y(-u) = for u < and Y(-u) = 
1 for u > 0. These expressions can be extended by continuity at u — ±0: 
we find X(a, b, +0) = 1 and Y(a, b, +0) = when u — > on the right, and 
= when w^Oon the left. Moreover, Eqs. (40)-(41) yield once again Eqs. 
(36)-(39) at u — ±1, provided that (gT)(a, u) is replaced by and the Cauchy 
integrals are replaced by ordinary integrals. 
Finally, putting z — > —1/k into Eqs. (34)-(35) leads to 



X(a,b,-l/k) = 1 + R(a, k) j[S(a, k) + kb][l - £ x (a, b, +l/k)\ 

+ ^ Y {a, b, +l/k) exp{-kb) + ^C( a > b > ~ l / k ) exp(+A;6)| 
+ § f J (M)j[l-f.W.)l I ^ + ««,W^j d„, (42) 

Y(a, 6, -1/Jfe) = -i?(a, fc) - £ x (a, 6, +1/Jfe)] exp(-A;6) 

b, -1/k) exp(kb) + [S(a, k) + kb]£ Y (a, b, +1/A;)} 

- 1 1" ^ «> {i 1 - ""•^i^rrlsr + « tt > "'r^} iv - (43) 

and putting z — > 

X(a, 6, +l/fc) = 1 + R(a, k) - £ x (a, 6, +1/A;)] 

+[S(a, fc) - kb]£ Y (a, b, +l/k) exp(-kb) - jg&(a, 6, exp(-A;fo) J 

+ | jf 1 <?M | [1 - & (a, 6, «)] + £ y (a, b, v f^'^ } dt;, (44) 



Y(a,b,+l/k) = -R(a, k) \[S(a,k) - kb\[l - £ x (a, b, +l/k)\ exp(-kb) 
+^(«. & > +V*0 exp(-A;6) + ^ y (a, b, +1/7V)} 
- | ^ <?M | [1 - & (a, M)] ^y^~ + Ma, 6, v) } dt;. (45) 
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It is possible to eliminate the derivatives at ±l/k of the functions £ x and 
£ y which appear in these four equations. Just use the derivative, at ±l//c, of 
Eqs. (16)-(17), remarking that T(a,±l/k) = 0, T'(a,±l/k) = ±k/R(a,k), 
£ y (a, b, — 1/k) exp(kb) = 1 — £ x (a, b, +l/k) and [1 — £ x (a, b, —1/k)] exp(kb) = 
£ Y (a,b, +l/k). The latter two identities follow from Eqs. (16)-(17) at —1/k. 
Therefore, at z — —1/k 

X(a, 6, -1/k) - (l/k)R(a, k)&(a, b, -1/k) exp(kb) = 

R(a, k) {kb[l - £ x (a, 6, +1/A;)] - (l/k)&(a, b, +l/k)} , (46) 

Y(a, b, -1/k) - {l/k)R(a, k)£ x (a, b, -1/k) exp(kb) = 

- R{a, k)[kb£ Y {a, b, +l/k)} + (l/kK' Y {a, b, +l/k)}, (47) 

and at z — +l/k 

X(a, b, +l/k) + (l/k)R(a, k)£ Y (a, b, +l/k) exp(-kb) = 

R(a, k)[—£ Y (a, b, +l/k)kbexp(—kb) + (l/k)£ x (a,b, -1/k)], (48) 

Y(a, b, +l/k) + {l/k)R{a, k)g x (a, b, +l/k) exp(-kb) = 

R{a, k) {[1 - £ x (a, b, +l/k)}kbexp{-kb) + {l/k)£ Y {a, b, -1/k)} . (49) 

Inserting these equations into Eqs. (42)-(45) leads to 

= 1 + R(a, k) js(a, k)[l - £ x (a, b, +l/k)) + ^' x (a, b, +l/k) 

+^ Y (a,b,+l/k)exp(-kb)} 

+ | jf 1 g(a, v) |[1 - £ x (a, 6, ^73^ + ^ (a ' & ' ^ T+fc^ } dw ' (50) 



= R(a, k) -^ x (a,6,+l/fc)]exp(-fc6) 

+S(a, fc)^(a, 6, +l/k) - jj-&(a, 6, +1 A)} 

+ £ jf <K«, «) {[1 - k(a, 6, «)l g YTto^ + ^ (a ' 6 ' ^T^} dv ' (51) 

= 1 + R(a, k) |i[l - £ x (a, 6, +1/*)] - ±& (a, 6, -1/k) 
+ 5(a, k)£ Y (a,b, +l/fc) exp(-fc6)| 
+ \ \l 9(a, v) |[1 - £ x (a, 6, + Ma, 6, ^L"^ } d ^ (52) 
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= R(a, k) | ^(a, k)[l-^ x (a,b,+l/k)]exp(-kb) 

+ ^ Y {a, b, +l/k) + i&(a, 6, -1/A;)} 
+ I jf 1 <y(a, v) |[1 - ^(a, 6, + ^ a > 6 ' ^T^jfeu} dw - (53) 

It is now possible to eliminate the derivatives appearing into Eqs. (42)-(45) 
by calculating them with the help of Eqs. (50)- (53), which gives 

1 • I k) - iAa.h. • 17./"' 



A;6 



X(a,6,-l/ife) = 1 + kbR(a, k) 

+ ^ (a ,M)exp(-6/„)l^i±™} d„, (54) 



y(a,6,-l/ife) = exp(A;6) 

+ fcM2(a, k) I [1 - ^(o, 6, +!/*)] - ^(a, 6, +1/*) 



A;6 



- 7^ g(a,v) |[1 -f x (a,&,u)]exp(-&/u)- — 



1 — exp[(l + kv)b/v] 



+ kv 



HAaM '-^-^ ldv, (55) 



X(a,6,+1/A;) = 1 

+ kbR(a, k) exp(-kb) |[1 - £ x (a, 6, +1/*)] - ^(a, 6, +1/A;)} 

a yi J l-exp[-(l + fc^)6/^] 

+ e y (a,M)exp(-^) 1 " eXP[ 1 '^ M ^ ] }d,, (56) 



y(a,6,+l/fc) = exp(-£;&) 
+ kbR(a, k) exp(-kb) 



l-S x (a,b,+l/k)-S Y (a,b,+l/k) 



sinh(A;6) 



I Jo ^ a ' u ^ _ ^( a ' fo '' f; )] ex P( _& / v 



A;6 

1 — exp[(l — kv)b/v] 
1 — kv 



/ , A - exp[-(l + kv)b/v]} , 
+£ y (a, 6, u) f dw - ( 57 ) 
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These expressions do satisfy relation (11), as it should. 

Some additional equations involving the functions £ x and £ y can be derived by 
letting z — > oo in the second equalities (34)-(35), using Eq. (22) and observing 
that T(a, oo) = 1 — a. These are 

1 - (a/2)a (a,b) = l+R ^ fc) ^ _ ^ ft> + ^ ft> exp( _ A;6) ] 



a 

+ 2 



/ ^(a,u)[l-^(a,6,u)+^(a,6,u)exp(-6/u)] du, (58) 
•/ o 



S " ''' 1 '° U ' - = R(a, k){[l- £ x (a, b, +l/k)} exp(-kb) + £ y (a, 6, +!/£;)} 



1 — a 



+ ^ f 1 g(a,v){[l-Z x (a,b,v)]exp(-b/v)+£ Y (a,b,v)} dv. (59) 
2 Jo 



Subtracting these equations from Eqs. (34)-(35) and their extension at ±1 and 
±l/h leads to a new series of integral equations for the functions X and Y 
(not written here). 

Equations in this series generalize the ones for the function if in a half-space, 
and Eqs. (27)-(31) are retrieved by letting b — > +oo into Eqs. (34), (38), (40), 
(50) and (58) respectively. Note that the if-function is nowhere involved in 
these equations. The auxiliary functions of finite spaces can thus be calculated 
without first deriving those for semi-infinite spaces. 

Equations (40)-(41) are the simplified version of two equations derived by 
Danielian [18] in application to Rogovtsov and Samson's calculation of the 
resolvent function in a finite slab [19]. These equations are, for u > 0, Eq. 
(40)[w] u>0 + exp(-&/w) x Eq. (41)[-u] u>0 and Eq. (41)[u] u>0 + exp(-6/u) x 
Eq. (40)[— u] u>0 , where Eq. (40)[u] u>0 denotes Eq. (40) with u > 0, and Eq. 
(40) [— u] u>0 denotes Eq. (40) with — u in place of u, supposing then u > 
[ditto with Eq. (41)]. As a result, Danielian's equations contain two times too 
many terms, a remark that still applies to its expressions of the functions / 
and / introduced for solving the transfer equation within a slab: see the end 
of the next section. 



5 Application to the albedo problem 



In this section, we consider the "restricted" albedo problem, which consists in 
finding the source function of the diffuse field within a slab illuminated from 
outside. The calculation of the specific intensity is left aside. It is supposed that 
parallel rays are incident on the top surface r = at an angle arccos uq to the 
inward normal {0 < Uq < 1) , with unit flux per unit area normal to the rays. 
Since the mean intensity of the direct field at level r is (l/47r) exp(— t/uq), 
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the source function for the diffuse field is given by 

S(a, b, t, m ) = T-B(a, b, r, u ), (60) 

where the function B = B(a, b, r, z) is solution to the following integral equa- 
tion 

B(a,b,T,z) = exp(-r/z) + ^ / £i(|r - t|)B(a, 6, t, z) dt, (61) 

2 Jo 

with parameter z G C*. i?i is the first exponential integral function as defined 
by E 1 (r) = Jo 1 exp(-r/«)( du/u) (r > 0). 

There are many analytical methods for solving the above equation. The oldest 
one, which is possibly the most direct, originates in Ambartsumian's solution 
to the problem of diffuse reflection of light in a half-space [8]. His approach 
was generalized by Busbridge to a finite slab [20]. A modern version of this 
method, which applies to complex values of the ^-parameter, is described in 
Sec. 3 of [21]. Ambartsumian's pioneer paper contains the first expression of 
the £?-function on the boundary layer of a half-space, viz. 

B(a,0,u o ) =H(a,u ). (62) 

In the finite case, it seems that the surface values of the _B-function have been 
given for the first time by Busbridge, who found in [20] 

B(a, b, 0, uq) = X(a, b, uo), B(a, b, b, uq) = Y(a, b, u ). (63) 

The 5-function may be calculated within the slab using the following relations, 
valid for any z G C* as seen in [21]. In a semi-infinite medium, Eq. (61) has a 
solution as long as z G C* satisfies the condition Re(l/z) > —k. This solution 
is 

B(a, r, z) = Y[- Re(z)] e *l[~ T/z) + H(a, z)zr)(a, r, z), (64) 

where the Heaviside function Y is extended by 1/2 at 0: Y(0) = 1/2. The 
77-function is defined for r > and z G C* by 

r}(a,r,z) = — H(a,l/z')exp(Tz') , (65) 

Z17T J-ioo 1 + Z Z 

and it is extended by continuity at r = [21]. 
In a finite slab, one has for z G C* 

B(a, b, t, z) = X(a, b, z)zr) x (a, b, r, z) — Y(a, b, z)zr) Y (a, b, r, z), (66) 

where 

Vx (a,b,r,z) = — / X(a,b,l/z')exp(rz')— —, (67) 

LYK J — ioo 1 + Z Z 

\ /■+ioo 

r] Y {a, b,r,z) = — / Y(a, b, l/z) exp(rz') - , (68) 

2l7l J-ioo 1 + z z 
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for any r g]0, b[ and z G C*. These functions can be extended by continuity 
at r = and r = b [21]. 

Integrals in the right-hand side of Eqs. (65), (67) and (68) are calculated 
along the imaginary axis in the sense of the principal value at infinity, i.e. 
f-ioo = limx^+oo I-ix ■ These are also Cauchy principal values when z is on 
the imaginary axis. All integrals can be transformed by the method of residues 
for the numerical evaluation of the functions r], i] x and r\ Y - In the semi- infinite 
case, details on this classical transformation can be found in [16] and [17]. 
In the finite case, the calculation is detailed in [17]. The results include the 
function g(a,v) and the coefficients R(a,k) and S(a,k) as defined by Eqs. 
(24)- (26). One obtains, in a semi-infinite medium 



. kR(a,k) expi—kr) AAFT ^ . exp(— r I z) 
r)(a, t, z) = — \ \ - K - — - + Y[Re(z)]H(a, -z)—^ — 



H(a,+l/k) l-kz 



+ U\g/H)(a,v)eM-rM— (69) 
2 Jo v — z 



forzGC\{[0,l]U{+l/A;}}, 



r)(a,T,+l) 



kR(a,k) exp(-kr) | a f 1 
H(a,+l/k) 



^ z ^ + ^j\ 9 /H)(a,v)eM-r/v)^- 1 (70) 



at z — +1, 



kR{a,k) exp(-kr) / v exp(-r/w) 

V{a,r,u) = — —77-— ; + {gT/H)(a,u)- 



H(a,+l/k) l-ku 



u 



+ U\g/H)(a,v)exp(-r/v)-^- (71) 

Zi J (J LI 



for z = u G]0, 1[, and 



77(0, r, +l/k) 



kR(a, k) 
H{a,+l/k) 



k H(a,+l/k) 

a rl 

~n k (9/H)(a,v)exp(-T/v) 
2 Jo 



exp(-kr) 
dv 



1 — kv 



(72) 



at z = +l/k. 

In a finite medium, we need only the expressions of the ?7 x -function, since 



r} Y (a, b, r, z) = i] x (a, 6, 6 - r, -z) 



(73) 
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due to Eqs. (67)- (68) and (11). These are: 
rj x (a,b,T,z) = 

kR(a, *){[!- Ua, b, +l/k)] e -^t^l - £ y (a, b, + 1/ k ^'f' " )] } 

+ {Ypte(z)] [1 - k(a, 6, *)] exp(-r/z) 

z 1 [a, z) 

+Y[- Re(z)]^(a, 6, -2) exp[(6 - r)/*]} 

+ ^ /'»(«■«) ([1 -U«.M1^^ - f,(a,M) ° Ph( ' '- T)M ) d» 
2 Jo [ i> — 2; f + zj 

(74) 

forzeC\{[-l,+l]U{±l/A;}}, 
rj x (a,b,r, -1) = 

fei2(a, fe) {[1 - k(a, 6, +l/k)]^t±± - £ y (a, 6, +l/fc)^!z^_Z^l j 

+ ;/^M)f[i-U a ,M)]^ + M fl ,M) exp| -f' r)H }d, 

2 JO [ 1 + 1> 1— fj 

(75) 



?7 x (a, 6,r,+l) = 

fci2(a, fc) { [1 - k(a, 6, +1/*)]^=^ - £ y (a, 6, +1 A) ^"^" ^ } 
- 1 f 9(a, v)[[l- Ua, b, v)] ^tlM + gy (fl> 6> w) «p[-(6 + - r)/t>] | ^ 

(76) 

at z — ±1, 
rj x (a,b,T,u) = 

kR(a, *){[!- £, («, 6, +1/*)]?S^ - MM, + l/t)S!z^Zl)l j 

+ i( 9 T)(o,«) {Y(u)[l -£*(<,, 6, u)]«p(-t/u) 

■u 

+Y(— u)£ y (a, b, -u) exp[(6 - r)/u]} 

+ l/'sC »)([!- Cv(a, 6, _ Wa , j, B) «P[-(t -r)/.] | 

2 Jo [ v — u v + u ) 

(77) 
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for z = u e] - 1, 0[U]0, +1[, and 




} 



dv, 



(78) 




r} x (a,b,T,+l/k) = -kR(a, k) \ [S(a, k) - kr][l - £ x (a, b, +l/k)] exp(-kr) 



at z = ±l/k. It is worth noting that these expressions only involve the func- 
tions C, x and £ y , not the functions if, X or Y. They can be checked putting 
in b — > +00, since then 7/ x does coincide with 7/ and r/ y vanishes. 
From the physical interpretation of the functions 77, ?7 X and 77^ we reminded 
in Sec. 5 of [21], it follows that these functions are continuous on the right 
at t = 0, and that r\ x and r\ Y are continuous on the left at r = b. Their 
restrictions to the boundary layers are thus given by the limits of the above 
expressions as r — > and r — > 6, which coincide with their values at r = and 
t = b. Using the integral equations of Sees. 3 and 4, one obtains the following 
surface expressions of the functions 77, 7/ x and i] Y , valid for any 2 in C*: in a 
semi-infinite medium 



dv 



(79) 



77(0, 0,z) 



1 



{l-Y[-Re(z)]if(a,-z)}, 



(80) 



and in a finite medium 



r) x (a,b,0,z) 



1 



{1 - Y[-Re(z)]X(a,6, -2)} 



(81) 



z 



r) x (a,b,b,z) = -Y[Re(z)]Y(a,b,z), 



(82) 



so that, from Eq. (73) 



77 y (a, b, 0, 2) = — Y[- Re(z)]Y(a, 6, -2), 



(83) 
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rj Y (a, b, b, z) = ~ {1 - Y[Re(z)]X(a, b,z)}. (84) 

These expressions still hold on the imaginary axis, since Y(0) = 1/2. Substi- 
tuting them into the relations (64) and (66) at r = and t — b, one obtains 
the following surface values of the 5-function: 

B(a, 0, z) = H(a, z) [Re(l/z) > —k], (85) 



B(a,b,0,z) =X(a,b,z) (zeC), (86) 

B(a, b, b, z) = Y(a, b, z) (z G C*), (87) 

which generalize in the complex plane the classical relations (62)- (63) derived 
over ]0, 1] more than half a century ago. The need for such an extension is ex- 
plained at the end of Sec. 2.3 of [22]. We conclude that the integral equations 
of Sees. 3 and 4 provide the connection between the internal and the surface 
solution to the albedo problem (61) with complex parameter z. 
In conclusion, we outline the bibliography on the exact solution of the albedo 
problem. We have seen that the surface solution (62) was given for the first 
time by Ambartsumian in a half-space [8]. The internal solution was first 
derived in the context of neutron transport by means of Case's singular eigen- 
function method [24,6]. It was expressed in terms of the r/-function by Danielian 
and Mnatsakanian [25] . The functions F and F of these authors are indeed re- 
lated to our function rj by F(a, r, u) = ur](a, r, u) and F(a, r, u) = ur](a, r, —u) 
(0 < u < +1). Similar functions have also been introduced by Crosbie and 
Linsenbardt [26] and Viik [27] . The first analytical expression of the ^-function 
is in Danielian and Pikichian [28], whose relation (36) is valid over [— 1, 0[, not 
over ]0, 1]. 

In the finite case, the surface solution (63) was given by Busbridge [20], and 
then Case's method was again required for calculting the internal field [29]. 
The first expression in terms of the rj x - and %-functions is due to Kagiwada 
and Kalaba [30,31], who used the invariant imbedding method. Their func- 
tions b and h coincide with our functions i] Y and r] x respectively. Danielian's 
approach in [32] is similar, since his functions F and F are related to r] x 
and i] Y by F(a,b,r,u) = uf] x {a,b,T,u) and F(a,b,r,u) = ur] x (a,b,r, —u) = 
ur) Y (a,b,b — t,u) (0 < u < 1). In [18], Danielian wrote the analytical expres- 
sions of his functions F and F, henceforth denoted by / and / in a finite slab. 
His formulae may be simplified with the help of Eqs. (40)-(41), since they 
contain two times too many terms. We recall that Danielian was not aware of 
Eqs. (40) and (41) separately (see the end of Sec. 4). After simplification, Eqs. 
(11)-(12) of [18] do yield our own Eq. (77) again. The latter and accompanying 
equations in the series (74)-(79) are given in [17]. 
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6 Conclusion 

The integral equations satisfied by the if-function (Sec. 3) and by the X- 
and F-functions (Sec. 4) are our main results. The former one generalizes, in 
the whole complex plane, an integral equation derived by Busbridge [1] over 
[— 1, +1]. The latter ones for the X- and y-functions generalize and simplify 
two integral equations derived by Danielian [18]. These integral equations 
clarify the slab albedo problem (61), since they show that its internal solution 
simplifies on the boundary layers into the surface solution. 
Contrary to Danielian [18], we consider that these integral equations are of 
little interest for computing the functions H, X, Y, £ x and £ Y . The if-function 
can be easily evaluated from one of its many analytical expressions. This 
approach is particularly interesting in strongly scattering media. The best 
way we know to compute the functions X, Y, £ x and £ y is that developed by 
Busbridge [1] and Mullikin et al. [33,34]. The four functions are derived from 
some simpler auxiliary functions that solve Fredholm integral equations with 
regular kernels. This approach is that of Danielian himself in two important 
papers [35,36] he published after Ref. [18]. It was also developed in Ref. [17], 
which contains ten-figure tables of the functions in question. 
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